Nonlinear patterns in Bose-Einstein condensates in dissipative optical lattices. 
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It is shown that the one-dimensional nonlinear Schrodinger equation with a dissipative periodic 
potential, nonlinear losses and linear pump allow for the existence of stable nonlinear Bloch states 
which are attractors. The model describes a Bose-Einstein condensate with inelastic two- and 
three-body interactions loaded in an optical lattice with losses due to inelastic interactions of the 
atoms with photons. 
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I. INTRODUCTION 

Being indispensable attribute of all natural processes 
the dissipation accompanied by properly adjusted com- 
pensation of losses, plays a constructive role in generation 
of nonlinear patterns, which nowadays arc observed in 
many natural sciences, including physics, biology, chem- 
istry, life sciences, etc [J. As attractors the respective 
structures play a prominent role in practical applications 
among which we particularly mention light patterns in 
nonlinear optics, receiving a rapidly increasing deal of 
attention [2| . For description of the emergence of dissi- 
pative patterns several mathematical models have been 
developed, the complex Ginzburg-Landau equation being 
the most simple and widely used one. This model allows 
for quantitative and often even qualitative description of 
a wide range of the physical phenomena [|[ . 

The dissipation is also an important factor in the the- 
ory of Bose-Einstein condensates (BECs), where it ap- 
pears in a natural way either through the inelastic inter- 
actions of light with atoms (see, e.g., Q), or in a form 
of inelastic two-body and three-body inter-atomic inter- 
actions [E IS 0- Alternatively, the dissipation can be 
introduced artificially, say in a form of nonlinear dissipa- 
tive lattice [1] or by probing the condensate by means of 
an electronic beam [9(, and can affect the matter wave 
dynamics in a variety of ways It is then natural to 
expect that a BEC where the dissipative losses are prop- 
erly compensated by a pump of atoms to the system, can 
give rise to highly stable atomic density distributions. 
Probably the simplest nontrivial and practically impor- 
tant realization of such situation can be a BEC loaded in 
a one-dimensional (ID) optical lattice (OL). 

This leads us to the aim of the present work: we report 
theoretical study of stable matter wave patterns emerg- 
ing in an array of BECs when the linear and nonlinear 
dissipations are compensated by the homogeneous pump 
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of the atoms. More specifically, we show that in such 
systems emergence of nonlinear dissipative Bloch waves 
is possible. Such periodic patterns appear as attractors, 
and hence initial distributions in a large range of the pa- 
rameters rapidly evolve to the dissipative Bloch waves. 
These are essentially nonlinear states, which in a general 
situation do not allow for transition to the limit of zero 
density, and at the same time, when the dissipation is 
not too strong, can be naturally associated with a spe- 
cific edge of the band of the underline conservative part 
of the lattice potential. 

The organization of the present paper is as follows. In 
Sec. QI] wc formulate the models of the dissipative cubic 
nonlinear Schrodinger (NLS) equation and provide the 
estimates for the main effects to be observable. In the 
next Sec. IIIII we perform numerical study of the emer- 
gence and stability of the periodic patterns. Sec. [IV] is 
devoted to extension of the results to the NLS equation 
with quintic dissipation. In Sec. [V] wc present the re- 
spective lattice models, which can be obtained using the 
tight-binding approximation from the initial continuous 
dissipative NLS equations. The results are summarized 
in Conclusion. 



II. THE MODEL AND ANALYTICAL 
ESTIMATES 

A. The model 

We start with the ID Gross-Pitaevskii (GP) equation 

i* t = -$> xx + 2asin 2 (a:)* + g\^\ 2 ^ + ir* (1) 

where a — a' — iot" is the properly normalized complex 
atomic polarizability whose imaginary part a" > be- 
comes appreciable for relatively large electric fields (see, 
e.g., [H), g = g' — ig" is the dissipative nonlinearity (with 
g" > describing inelastic interatomic collisions) and 
r > is the linear gain. 

It is relevant to mention, that emergence of coherent 
localized and periodic structures from the interplay be- 
tween periodicity and complex (i.e. conservative and 
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dissipativc) nonlincarity have recently being addressed 
within the framework of the complex Ginzburg-Landau 
equation [ill ] and nonlinear dissipativc Bloch waves in 
a model of optical parametric oscillators were reported 
in Ref. fl2| . Compared to these previous studies, apart 
from the new physical applications, the present model 
(U) display very different mathematical features, which 
stem from the periodically varying dissipation which is 
controlled by the OL. 

For the next consideration it is convenient to introduce 
the parameter 



a 



(2) 



which will play the role of the effective gain (see below), 
and parameter <f> defined by the relations g"/\g\ = sincf) 
and g'/\g\ = cos</), which controls the nonlinear dis- 
sipation (the parametrization introduced implies that 
< <t> < 7r). Then, renormalizing the macroscopic wave- 
function -0 = \g\^l 2 e~ %a we rewrite Eq. ([T]) in the form 



ty t = -^ xx -(a'-ia")cos(2x)i; + e'" l4 '\ip\ 2 ^ + iA^. (3) 



B. On homogeneous linear pump 

The model (J3]) explored in the present paper is based on 
the supposition that there exists some " gain" mechanism 
compensating dissipativc losses. Generally speaking such 
a mechanism would require an infinite external reservoir 
of atoms supplying the condensate. Inclusion of the linear 
gain has been explored in the BEC applications [H, [l3| 
and is usually associated to the growth of the number of 
condensed atoms due to condensation from the thermal 
cloud [f|. However, a possibility of controlled and a priori 
given gain of the number of atoms, to the best of our 
knowledge, has not been proposed, so far. 

In this subsection we outline an idea on how an ef- 
fective homogeneous gain can be implemented " mechani- 
cally" . It is based on the fact that, what really matters is 
the linear atomic density rather than the total number of 
atoms. Hence, one can consider a cigar-shaped conden- 
sate whose dimensions decrease in time, thus resulting in 
increase of the linear density. 

To this end we consider a ID BEC loaded into a 
cigar-shaped potential, which consists of a periodic OL 
and a parabolic trap. In order to implement an effec- 
tive linear homogeneous pump we assume that (i) both 
the longitudinal and transverse harmonic oscillator fre- 
quency are growing functions of time: wm(T) = u>i\f(T) 
and luj_(T) = ujj_f(T), where u>u and uj± are the re- 
spective frequencies at initial time T = and f(T) 
is a positive definite increasing function, chosen to be 
f(T) = (1 — 2(jj||T)~ 1 ; (ii) the complex s-wave scatter- 
ing length grows in time: a(T) = a s f(T), a s = a' s — ia" 
being the initial complex scattering length (the desired 
time dependence can be achieved with help of Fcsh- 
bach resonance), (iii) the lattice amplitude and the lat- 
tice period are functions of time: V(T) = Vf(T) and 



d(T) = df-^ 2 (T), V = V - iV" and d being the re- 
spective initial values (this can be achieved by changing 
the laser beam intensity and the angle between the beams 
creating the OL). Then, in the mean- field approximation 
such a condensate is described by the 3D GP equation 
with varying coefficients: 



iM T = -— V 2 * + 2V(T) sin : 



in 

T 



2m 

-2 , , ,2 



Gen) 



to 2 (T)X 2 + uj 2 l (T)R 2 ± 



4irh 2 a(T) 



*| 2 *, (4) 



where ^ and m are wavefunction (in the physical units) 
and the mass of the atoms, respectively, = (Y, Z) is 
the transverse coordinate. 

In the case when OL period d is much bigger than 
the transverse harmonic oscillator length, d 3> ftj_ = 
\Jhj (tocjjJ, Eq. ([4]) can be reduced to the ID GP equa- 
tion by means of the multiple-scale expansion proce- 
dure [lj]. To this end the 3D bosonic wavefunction has 
to be searched in the form 



*(X,R X ,T) = 



^#M)C(R,,T)x 



exp | -i—^- - i— log[/(T)] - i \ _ ) , (5) 



2uj« 



2u±a 2 l 



where ip(x,t) is the unknown dimensionless function of 
the new dimensionless independent variables 



= -v//(T)X 



and t = 



2hcu, 



log[/(T)] (6) 



(hereafter E r = ti 2 n 2 / (2md 2 ) is the "initial" recoil en- 
ergy) and C(R±,T) = 7T- 1 / 2 a- ± 1 cxp[-Rlf(T)/(2a 2 L )] 
describes the linear transverse distribution varying in 
time due to change of the transverse trap. Then, in- 
troducing the dimensionless parameters a = a' — ia" = 
{V - iV")/E r , A = 3hw ll /{2E r ) - a", and 0, the lat- 



ter defined by the relations a' s ' 



a, = sin i 



and ai 



cos 4>, and substituting the ansatz © in Eq. (H} one read- 
ily obtains the desirable equation 

As it is clear the above mechanism formally works only 
for T < l/2w|| and implies change of the trap dimen- 
sions, which should not affect the applicability of the 
meanfield model (HI- This means that the created ef- 
fective pump is relatively small, and thus can properly 
work for small dissipative losses. Say, if in the dimension- 
less units r ~ a" ~ 0.001 , for the rubidium condensate 
loaded in a cigar-shaped trap with the initial dimensions 
a±_ = 0.5 fim and a\\ = 17.32 /im (corresponding to the 
u)h w 2.42 Hz), and with the imposed optical lattice with 
the constant d = 1 /zm, one computes that the time in- 
terval during which the constant pump can be supported 
is T « 155 ms (or t ft 1040 in dimensionless units), and 
corresponds to the decrease of the linear dimensions of 
the trap two times. 

Below, in this work, however we explore more strong 
values of the gain (and dissipation), in order to reduce 
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time during which instabilities are developed and unsta- 
ble initial conditions converge to the respective attrac- 
tors. More specifically, following [(| for all numerical 
simulation we will chose T = 0.1. 



C. On the spectrum of the linear problem. 

We are particularly interested in periodic stationary 
patterns. Therefore we explore the ansatz ip(x, t) = 
<p(x) exp(— i[it), with a real chemical potential (i leading 
the nonlinear eigenvalue problem 



flip = L<p + i[A + a" cos(2x)]ip + e~ i(t> \p\ 2 <p 



(7) 



,2 

where we have introduced the linear operator L = — — 
a'cos(2x). The eigenvalues and eigenfunctions of this 
operator are obtained from the Mathieu equation 



Lifink = E n {k)tp nk . 



(8) 



It gives the band-gap spectrum, characterized by the 
number of the band n > 1 (n = 1 corresponding to 
the lowest band) and by the wave vector k in the first 
Brillouin zone, k G [—1,1]. Subsequently, the spectrum 

of the operator L is given by E n (k) <G [e[~',e[ + '] U 



(E^ < E { + ] < E { - ] < 



Here 



[Et\4 +) ]U- 
En "* and E n + ^ are respectively the lower and the upper 
edges of the ra-th band, E n ~ ] = E n (0), E n +) = E n (l), if 
n is odd, and Ei~' = E„(l), E n + ^ = E n (0), if n is even. 

The requirement for \i to be real readily gives the re- 
lation 



[a" cos(2cc)|v?| 2 - sin(0)|^| 4 ] dx + AN = (9) 



where N = \i(;\ 2 dx is the normalized number of atoms 
per one lattice period. We notice, that the condition © 
can also be obtained from the requirement for the number 
of condensed atoms to conserve, i.e. from dN/dt = 0. 

Our analysis will be restricted to the case of weak dissi- 
pation, which we express in terms of the small parameter 
e = \J\a" ja'\ -C 1 and respectively we require A = e 2 S 
with \S\ < 1. It is then natural to recall that in the 
absence of the dissipation Eq. © possesses branches 
of periodic solutions which in the linear limit, i.e. when 
N — > 0, are reduced to the conventional Bloch states <p n k- 
Moreover, when N <C 1 the modulational instability of 
the respective periodic solutions is described in terms of 
the multiple-scale approximation ■ Therefore we start 
the description of the dissipative problem at hand with 
this limit. 

As the first step we consider the case when in the re- 
spective linear problem 



Lip + ie 2 V(x)ip = Eip, 



(10) 



the dissipation V(x) = 5 + a' cos(2a;) is treated as a per- 
turbation. It is straightforward application of the per- 



turbation theory to obtain corrections for a chosen eigen- 
value E n {k) due to the dissipative term: 



with 



and 



E = E n (k) + e 2 E^+e 4 E^ 



E { nl=i V{x)\ Vnk {x)\ 2 dx 



(11) 



(12) 



E 



(2) 
n k 



\Io V{x)pnk(x)p n , k (x)dx\ 2 



E n {k)-E n .{k) 



(hereafter an overbar stands for the complex conjuga- 
tion). The obtained result has several important conse- 
quences. Indeed, for the eigenvalue to be real, one has to 
require E^ = 0, or explicitly 



A + a"-f = 0, 7 



cos(2x)\<p nk (x)\ 2 dx. (14) 



This requirement for the eigenvalues to be real either can 
not be satisfied or can be satisfied only for a single gap 
edge [since 7 depends on the indexes (n, k) which are 
omitted for the sake of brevity of notations] . Thus in a 
general situation A / 0. If A < —a" 7, then E^) de- 
scribes effective dissipation, what means that the zero 
solution (if,' = 0) is an attractor. If however A > — a"j, 
the zero solution becomes unstable and one should ex- 
pect emergence of nonlinear coherent structures, which 
of course are possible if the respective linear pump is 
compensated by the nonlinear dissipation. This is pre- 
cisely the case of sin (f> > in the paramctrization chosen 
above [see also Eq.©]. 

For A = (this is the situation explored in all nu- 
merical simulations reported in this paper) the above ar- 
guments mean that the nonlinear stable structures are 
expectable only for 7 > 0. Then nonlinear patterns are 
generated form arbitrarily small initial fluctuations (see 
also Fig.[2]below). Meantime the sign of 7 is not uniquely 
defined and depends on the band-gap edge. For the cos- 
like lattice, explored here, positive 7 is verified only for 
the lowest bands. Hence, only close to the lowest band 
edges one can expect to obtain low-density stable nonlin- 
ear patterns. This conjecture is verified for all numerical 
simulations performed in the present paper. 

Finally, we observe the validity of the inequality 

| 7 | < 1 (15) 
which will be used for the discussion in what follows. 



D. Multiple-scale analysis of the stability 

By varying the pump on can selectively excite periodic 
patterns, whenever more than one stable state exist. In 
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the low-density limit, i.e. at N — > 0, the stability of 
the periodic solutions can be tested using the multiple- 
scale analysis, similarly to the approach successfully used 
in the conservative case Now e has to be used as a 
small parameter of the problem. Following this approach, 
for the chosen band edge, say E n , the wavefunction 
is approximated by ip w £ A({;, r)ipn (x) exp (— iE n t) 
where t) is an amplitude depending on slow vari- 
ables £ = ex, r = e 2 t, and ip^^x) is the Bloch wave- 
function at the edge of the n-th band. We observe, that 
by choosing e to be the small parameter we implicitly im- 
pose the conditions where the characteristic scale of the 
excitations is determined by the complex part of the po- 
tential (contrary to case of a conservative system, where 
the small parameter of expansion is treated as a free pa- 
rameter related to the detuning of the chemical potential 
towards the adjacent gap, see e.g. flij). 

Using the standard algebra (the details can be found 
say in Ref. 14[) one verifies that A(£, r) solves the dissi- 
pative nonlinear Schrbdinger equation 

iA T = -(2M)- 1 A^ + iAA + e^xl^A. (16) 
Here M = (d 2 E^ /dk 2 )' 1 is the effective mass, 



A = a'7 + S = 



a"7 + A 



(17) 



is the effective gain (according to the analysis of the 
Sec. Ill CI nonlinear patterns emerge near edges with 
A > 0, the conclusion which is also confirmed by the 
next analytical consideration and numerical simulations, 
what justifies the definition of A as a gain), 



X 



\<p 



^\ A dx 



(18) 



is the effective nonlincarity, and 7 is defined in (|14|) with 

<Pnk = </'n ± ' ) • Within the framework of the multiple-scale 
expansion the condition ([9]) for conservation of the num- 
ber of atoms can be recast as follows 



A - sin(0) X |A| 2 =0. 



(19) 



Eq. (|16p possesses a stationary, coordinate indepen- 
dent, solution in the form 



-4, 



A 



Xsin(0) 



(20) 



Here f2 = A/ tan(0) is a constant determining the shift of 
the chemical potential outwards the gap edge: recall that 
now /1 = E^ ±} + e 2 Q. Unlike in the case of conservative 
systems, the value of fl is fixed. It is determined by the 
linear density of atoms and by the nonlinear dissipation 
controlled by the parameter <f> through the relations 



N 



a" j + A 
X sin 4> 



M = 4 



a"7 + A 
tan(0) 



(21) 



Several important conclusions follow immediately from 
Eqs. ([21]) . First, the limit N — > can be reached only 
when the condition (fT4|) is satisfied, i.e. only subject 
to the requirement for the spectrum to be real. In this 
limit /.t — > E n . For the general case one must have 
N > what leads us to the condition on the linear 
dissipation a""f + A > necessary for existence of the 
nonlinear periodic solutions. This constrain, which in 
Sec. Ill CI was obtained form the analysis of the linear 
spectrum, is equivalent to the requirement A > in 
Eq. (|16p , and under the chosen parametrization of nonlin- 
earity phase cf> is necessary for the particle conservation 
(fT^j) . In the original notations, the obtained necessary 
condition reads T > (1 — j)a" and means that the pump 
should be stronger than the dissipation a" reduced by 
the lattice factor (1 — 7) [notice that due to (fl~5f this 
is a positive factor]. Now the number of particles per 
one OL period cannot be smaller than the minimal value 
N m in = (ce""f + A)/x, which is the number of particles 
in the case of pure dissipative nonlinearity, i.e. when 
(0 = tt/2). In this limit fi -> £& ±5 . 

As the second relevant property of the system as hand, 
following from (f2Tj) we emphasize that the chemical po- 
tential is not a free parameter, as this would happen in 
a conservative system, but is determined by the balance 
of incoming and dissipating atoms. This balance can be 
controlled say by the nonlinear dissipation, parametrized 
by 4>. Respectively, obtaining solutions corresponding to 
different branches (see e.g. Fig.Q] panels b and d) require 
change of the nonlinear dissipation. In practical terms 
this can achieved, say, by using variations of the light 
with a frequency close to a resonance with one of the ex- 
cited atomic state. Then the control parameter between 
the real and imaginary parts of the scattering length is 
the relation between the frequency detuning 8 form the 
excited level and the natural line width T: S/T [6j, and 
practically any value of the parameter (f> is achievable 
(see, e.g., the examples considered in 0). 

As it is customary, to check the stability of 

constant amplitude solutions, we consider small 

perturbations of the stationary solution: A st + 
^ ae iki-iwr + -^ e -iki+i^ g-inr ( here | a |j 6 | < 

Substituting this ansatz in Eq. (|16p and linearizing with 
respect to a and 6, we obtain a dispersion relation in the 
long- wavelength limit (k 1) 



-iA ± yJVlk 2 /M - A 2 



(22) 



Thus for the stability of the background A we have to 
require two conditions: 



A>0 and Ou - E&)M > 0. 



(23) 



The first of these inequalities has simple meaning of the 
positive effective dissipation which must compensate the 
nonlinear losses. The second condition is of " conservative 
origin" and coincides with the stability of Bloch states in 
the conservative systems (see, e.g., |14I|V 



5 



Besides the stationary solution (|20|) . Eq. ([T6]) possesses 
a time-dependent homogeneous solution 

A(t) = a(r) cxp (-i cos((f>)x J a 2 (T')dr'^j, (24) 



where 
o(r) = 



a 



^ Ae 2A(T ~r) + xs i n (^) a 2 (1 _ e 2A(r -r)) 



ao = a( T o)i and To is the initial time. 

For a positive A and nonzero ao the distribution (|24|) 
at r — > oo converges to the stationary solution (|20|) . 
In other words, the pattern ([20]) is an attractor for 
any nontrivial initial distributions, corresponding to a 
smoothly modulated Bloch function. The respective 
time-dependence of the particle number per period can 
be represented as 



TV 



N (a"-/ + A) 



N xsm(j) + e 2 ( Q "7+A)(to-i) (q/' 7 + A - iV xsii 



(25) 

where Nq is the initial number of particles, and we re- 
stored to the original dimensionless variables. 

Below we will also explore the dynamical regimes when 
after some interval of time the pump is switched off, i.e. 
r = and A = —a" (respectively, 5 = —a'). Recalling 
the definition (fT7|) as well as the property (fT5*|) . we con- 
clude that this is the case of A = a' (7 — 1) < 0. Now 
the particle conservation condition (fT9|) can not be met 
for the dissipativc nonlincarity sin(</>) > 0. It is evident, 
that in this case nonstationary solution ([24)) evolves to- 
wards zero, i.e. the zero solution becomes an attractor. 
The respective dependence of the number of particles per 
period on time is now expressed as 



N(t) 



Ar a"(l- 7 )e 2Q "(T- 1 )(*- tt ») 



a"(l - 7) + AT o Xsin(0) (l - e 2«"(7-i)(t-*o)) 



(26) 

At t — to ^> [a"} this decay has the exponential asymp- 
totic 



N(t) ~ iVott"(l-7) 2a "(7-i)(t-t ) 
W «"(l-7) + Xsin(0)7V o 



(27) 



III. NUMERICAL STUDY OF COHERENT 
STRUCTURES 



When the external parameters of the system (i.e. a', 
a", A and 4>) are fixed, with each band-edge one can 
associate only one solution, having the number of par- 
ticles and the chemical potential determined from (|2~lj) 
and varying when either of the parameters is changed. 
Bearing in mind physical applications of the model we 
concentrate on the study of the dependence of the num- 
ber of particles on the chemical potential N(fi) [what is 



the same as N(fl) due to the link between f2 and /z]. The 
respective analysis was performed for the lowest bands 
and is summarized in Fig. [TJ 

The first lowest band (n = 1) is associated to two 
branches of periodic solutions of Eq. ([3]) depicted in FigUJ 
(panels a and b). One branch is 27r-periodic (the branch 
A-B) having zeros of the density in maxima of the OL in 
Eq. |JJ), i.e. x p = 7r/2 + pit (where p is an integer). For 
this branch the dependence N(fi) reaches its minimum 
in the vicinity of e[ + \ The second branch is 7r-periodic 
(the branch C-D). It has no zeros in the whole space, 
this reflecting the fact the respective linear Bloch wave 
of the underline conservative system has no zeros. The 
dependence iV(/i) of the second branch reaches minimum 
at E[ \ We notice, that analytical expressions for N(p) 
and 4>(ia) in the vicinity of remarkably well match 
the numerical ones (respectively dashed and solid lines 
in Fig.EEKb). 

Turning to the upper edges, in Fig. [TJ;,d we show the 
branch E-F, which is associated with the upper edge of 
the second band. It has the minimum of the number of 
particles in the vicinity of E^ given by N m i n ps 0.026, 
i.e. smaller than the one in the vicinity of first band edges 
{N m in ~ 0.1, see FigQJi). We also clearly observe that 
the minimum of iV(/z) is considerably shifted to the lower 
values of (i in comparison with the analytical prediction 

(EID. 

Passing to the stability of the solutions, direct numer- 
ical integration of Eq. ([3]) reveals that 27r-periodic pat- 
terns are stable at fx < e[ + ^ (solution A in FigUjt,b) and 
unstable at fi > E^ (solution B in Fig. [TJi,b) . As this 
is typical for dissipative systems, the stable solution is 
an attractor, and therefore an unstable initial condition 
rapidly evolves towards the stable one having the same 
nonlinear phase 4>. This process is illustrated in the panel 
B of Fig[TJ where the unstable solution transforms into 
the attractor corresponding to the point D (c.f. panels A 
and B in Fig. |T|) . At the same time 7r-periodic solutions 

are stable at p > E[~^ (see, e.g., solution D in Fig0Ji,b 
and respective panel D) and unstable at fj, < e[ ' (sec 
e.g., solution C in Fig.[T^,,b and respective panel C, where 
we show the convergence of the unstable solution to the 
attractor corresponding to the point A). 

The described stability properties well corroborate 
with the simple analysis performed in the preceding sec- 
tions. This is not the case, however, for the solutions 
associated with the upper band edges. They appear to 
be dynamically unstable, again rapidly converging to the 
stable lower branches having the same parameter <p. This 
is illustrated by the evolution of the patterns correspond- 
ing to the points E and F shown in the respective panels 
of Fig. [TJ where they transform into the solutions, corre- 
sponding to points A and D, respectively. This however 
should not be viewed as the discrepancy with the above 
analysis of the stability. Indeed the latter one is done 
for slowly varying perturbations thus representing only 
necessary but not yet sufficient condition for the stabil- 
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FIG. 1: (Color online) The number of particles N (panels 
a,c) and the phase of the nonlinearity, <f>, (panels b,d) vs the 
chemical potential fj, for a 1 = 3.0, a." = 0.1, A = 0.0 obtained 
analytically (dashed lines) from Eqs. (|21[) , and numerically 
(solid lines). Grey strips represent the first (in panels a,b) and 
the second (in panels c,d) bands. The lower panels A-F illus- 
trate temporal evolution of periodic solutions corresponding 
to the respectively labeled points of upper panels. The points 
A, C, and E are chosen to have the same nonlinearity phase 
tj> « 3.055. Similarly, the points B,D, and F are characterized 
by the phase <f> « 0.116. The corresponding parameters 7 are 
7 = 0.641 (panels A,B), 7 = 0.536 (panels C,D) , 7 = 0.123 
(panels E,F). 



ity. Moreover the reported behavior well agrees with the 
conclusion about emergence of complex eigenvalues due 
to weak dissipation [see (fl"4]) as well as the subsequent 
discussion] . 

Having discussed the stability properties of the period- 
ical patterns we address the two natural questions: how 
one can generate the above-mentioned stable periodical 
patterns, and how they evolve if after some time the gain 
is off? The answer to the first question is rather simple, 
since the final pattern is an attractor it will be generated 




FIG. 2: (Color online) (a,b) The evolution of the density, 
obtained from the numerical integration of Eq. (j3} with the 
initial condition ip(x, 0) = 0.01. (c,d) The evolution of the 
density started with the stable patterns corresponding to the 
points A and D of Fig. [1] with the gain turned off at to = 500. 
(e,f) Evolution of the number of particles per one period cor- 
responding to the dynamics shown in panels c,d, respectively, 
calculated numerically (solid curves) and analytically from 
Eq. (|26[) (dashed curves). Notice, that solid and dashed lines 
are indistinguishable in the scale of panels e,f. The other pa- 
rameters are a' = 3.0, a" = 0.1, A = 0.0, <j> ~ 3.055 (panels 
a,c,e) and (f> * 0.116 (panels b,d,f). 



from any initial distribution having small density. This 
is illustrated in Fig. [5^,b where one observes evolution 
of a homogeneous nonzero excitation towards the stable 
periodic pattern, whose parameters (number of particles 
per OL period and chemical potential) are determined by 
the fixed phase <fi of the nonlinearity. More specifically, 
in Fig. the initial excitation is transformed into the 
solution, corresponding to point A in Fig. [T] while the 
initial excitation in Fig. evolves towards the solution 
with the parameters of the point D in Fig. [TJ It is worth 
to notice, that such a spatio-temporal evolution confirms 
the behavior, qualitatively predicted by Eq. (|24| . 

The answer to the second of the above questions be- 
comes evident from Fig. [2}>f, where wc show the evo- 
lution of the stable periodic patterns, after the gain is 
turned off (recall that turning off the gain is equivalent 
to setting A = —a"). Thus, after turning off the gain, the 
condensate starts losing the atoms and its density homo- 
geneously goes to zero [as it was predicted by Eq. 
The panels e and f of Fig. [2] demonstrate remarkably 
good correspondence between the numerical results and 
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the analytical predictions. 

IV. COHERENT STRUCTURES WITH 
QUINTIC DISSIPATION 

Significant dissipative losses in a BEC can occur due to 
the inelastic three-body interactions @, @, 0] ■ Therefore 
as the next step we address their effect on emergence of 
the coherent structures in the dissipative model 

i$! t = -& xx + 2asin 2 (a;)^' + g'\^f\ 2 ^f + ir* 

-i.g"|*| 4 *, (28) 

where g" > 0, and other parameters are defined as above. 
Similarly to the previous case, introducing the parame- 
ters A = r — a", g = g'/^/g", and the renormalized 

macroscopic wave-function ip = (g") 1 e~ ia t \f', we ar- 
rive at the equation 

ii\) t = —ip xx — (a' — ia") cos(2x)ip + iAip 

+.g|^| 2 ^-z|V|V- (29) 

Now the conservation of the number of atoms requires 
[c.f. Eq. ©] 

/■7T 

/ [a" cos(2x)\ip\ 2 - \ip\ 6 ] dx + AN = 0. (30) 
Jo 

Following the approach developed above for the 
analysis of the cubic dissipative term, now we 
employ the multiple-scale analysis of Eq. (j2"5)) 
representing the wavefunction in the form 
ip w e 1 ^ 2 A(£ : ,T)(p ± (x)cxp(~iEi ±) t). The slowly 
varying amplitude A now is governed by the equation 

iA T = -(2M)" 1 % + iAA + G X \A\ 2 A - iT\A\ A A, (31) 

where T = J Q |(/3 ± | 6 dx, g = eG (T ~ \G\ ~ 1) and other 
parameters are defined as in (|16[) . We emphasize that 
unlike in the cubic case, here we imposed the condition 
of smallness of the conservative nonlinearity coefficient g, 
since in this case the effects of the elastic two-body inter- 
actions and the inelastic three-body interactions become 
of the same order and must be accounted simultaneously. 
Now the constant-amplitude solution reads 

A st = (A/T) 1/4 e- 4 ° T . (32) 

The number of particles per OL period N and the chem- 
ical potential \x depend on the nonlinearity g 

j a" j + A (±) / a" j + A 
= V T ' M = ™ 5X V T ' ( ' 

As it can be seen from Eqs. (|33| . in the case of dissipa- 
tive quintic nonlinearity the analytically estimated den- 
sity of particles in a stationary mode does not depend 
upon the cubic nonlinearity g [unlike this happens in the 



cubic case, see (|2~Tj) ]. The two-body interactions, how- 
ever, determine the chemical potential /x for given N. 
In particular, since x > [see the definition (|18p ] we 
have that /i < (fj, > E^) when g < (g > 0). 

Fig. [3] illustrates that the simple analytical estimate for 
N(fj) (shown in panels a,c by dashed lines) is in agree- 
ment with the numerical values only in the vicinity of 
the band edges (unlike this was in the case of the cubic 
dissipative nonlinearity where the domain of quantitative 
coincidence of the analytical and numerical results was 
significantly larger, see Fig. [1]). The analytical estimates 
for g(fi) is in much better agreement with the dependence 
found numerically (Fig. [3j3,d). 

Turning to the simple stability analysis of the dissipa- 
tive Bloch waves, first we notice that the analysis per- 
formed in Sec. Ill Cl continues to be valid. Moreover, ana- 
lyzing the stationary solutions (|32|) we arrive at the con- 
ditions (|2"3"|) , earlier obtained for the cubic case. Also sim- 
ilarly to the cubic case, direct numerical simulations con- 
firm the validity of (|2"3"|) for the first band (see Fig. [3^i,b, 
where modes A and D are stable while modes B and C are 
unstable) and their failure for the top of the second band 
(see Fig.|3J:,d showing that modes E and F are unstable). 

V. ASSOCIATED DISSIPATIVE LATTICES 

As the last point, we observe, that the limit of a rela- 
tively large potential depth can be described by the tight- 
binding approximation which is obtained by discrctizing 
Eq. using the expansion 

lj)(x,t) = 22 c nk(t)Wnk(x) (34) 
1 1 k 

over the basis of the Wannier functions 

w nm (x) = -L / <p nq {x)e-™ m <>d q . (35) 
V2 J-oo 

For more precise conditions of the validity of the approx- 
imation as well as for the details of the expansion we re- 
fer to [IH , only mentioning here that it must be required 
sufficiently fast convergence of the coefficients uj nm of the 
Fourier expansion of the energy: 

i r 1 

w„o > w n i > ' ' ' , u mn = -J E n (q)e " mq dq (36) 

allowing one to restrict the consideration to the hopping 
of only nearest neighbors. 

At this point however we emphasize one important fea- 
ture which distinguishes the Wannier function mapping 
of the continuous model to a discrete one in the dissipa- 
tive case. While, in general, the tight-binding limit does 
not account a number if important features of the dy- 
namics of the underline continuum model [l5[, the tight 
binding approximation can be formally performed for a 
number of the lowest bands (provided the depth of the 



8 



0.24 
0.22 



c 

E 


F 






3.6 


3. 




4 


d 
















3.6 3. 




4 




27t o 



2it o 



FIG. 3: (Color online) Number of particles N (panels a,c) and 
cubic nonlinearity g (panels b,d) vs chemical potential /j, for 
a' — 3.0, a' — 0.1, A = 0.0 calculated analytically (dashed 
lines) from Egs. p3[) , or numerically (solid lines). Grey strips 
represent the first (in panels a,b) and the second (in panels 
c,d) bands. In the lower panels A-F we show temporal evolu- 
tion of periodic solutions at respective points of upper panels. 



potential is high enough), still giving reasonable approx- 
imation for the static solutions. In the dissipative case 
the situation is dramatically changed because of the sta- 
bility properties. Namely, now the tight-binding approx- 
imation makes sense only for a band where the periodic 
solutions are attractors. Even when other bands result 
in discrete dissipative models with stable solutions, dif- 
ference between such solutions and those of the original 
continuum model grows exponentially already at early 
stages of the evolution. Therefore, in what follows we as- 
sume that the expansion (|34[) is performed for the band 
associated to the stable solutions. Moreover, in accor- 
dance with the examples considered above we choose the 
first lowest band with n = 1. This assumption allows us 
to drop the band index n below. 

Then, using the orthogonality of the Wannier functions 



we readily obtain the dissipative discrete NLS equation 
for the exnansion coefficients c„ 



for the expansion coefficients 
( a"\ , 

I 1 ~ i_ ^7 I (^oc m +wic m+ i +u>ic m _i) 

+e-^W\c m \ 2 c m . (37) 



dr 
! — — 

dt 



wf m (x)dx and we have neglected hop- 



Here W = f 

ping integrals with the upper bands and among next 
nearest neighbors. The obtained equation (j3~7| was con- 
sidered earlier in (l6j (it is also relevant to mention stud- 
ies of a more general discrete dissipative model in [171]). 
Therefore we do not proceed with the further analysis of 
(f3"T)) . referring for the study of their modulational insta- 
bility to the mentioned work. 

We complete this last section with the indication that 
the qubic-quintic dissipative model (|29p can be dis- 
cretized in the similar manner 



. dc„, 
' dt 



( «"\ , 

= ( 1 - 4-^7 I (UJoCrn + UlC m +l + LOiC m -l) 



+gW\c m \ 2 c m -iB\c m \ i c m , (38) 



where B =_J^' 00 w\ m {x)dx. This model was also ad- 
dressed in uM. 



VI. CONCLUSIONS 

To conclude, we have reported emergence of nonlinear 
periodic structures in a nonlinear Schrodinger equation 
with a dissipative periodic potential, nonlinear losses and 
linear pump. The models describe a Bosc-Einstein con- 
densate loaded in an optical lattice, where the potential 
created by the light interaction with the condensate ac- 
counts for nonelastic interactions of photons with atoms. 
We addressed nonlinear losses due to two- and three- 
body interactions. These simple one-dimensional model 
reveals existence of stable nonlinear Bloch waves (attrac- 
tors) to which large range of the initial data converges. 
The obtained solutions do not have the linear limit (i.e. 
the limit of the zero density) and the number of particles 
is limited from below. In all the simulations performed 
the observed stable periodic patterns correspond to the 
lowest band. The existence of the attractor has impor- 
tant physical consequence: in the described models load- 
ing particles (using the linear pump) is only possible to 
the given state, which is imposed by the nonlinear losses. 
By changing the interatomic interactions one can follow 
different branches of the stable solutions. 

Meantime the simplest stability analysis performed in 
the present paper does not yet describe all the features of 
the system which are related, in particular, to the scales 
comparable with the lattice constant, and thus not ac- 
counted by the multiple-scale expansion. These issues of 
the theory are left for further studies. 
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